arXiv: 1501.0372lv2 [cond-mat.supr-con] 28 Mar 2015 
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We report on highly tunable radio frequency (rf) characteristics of a low-loss and compact three 
dimensional (3D) metamaterial made of superconducting thin film spiral resonators. The rf trans¬ 
mission spectrum of a single element of the metamaterial shows a fundamental resonance peak at 
~24.95 MHz that shifts to a 25% smaller frequency and becomes degenerate when a 3D array of such 
elements is created. The metamaterial shows an in-situ tunable narrow frequency band in which the 
real part of the effective permeability is negative over a wide range of temperature, which reverts to 
gradually near-zero and positive values as the superconducting critical temperature is approached. 
This metamaterial can be used for increasing power transfer efficiency and tunability of electrically 
small rf-antennas. 


INTRODUCTION 

A three dimensional (3D) array of sub-wavelength 
elements composing a metamaterial can demonstrate 
unique electromagnetic properties which are not directly 
accessible in nature, such as negative permeability mm 
and zero refractive index 0 E]. Most metamaterials 
achieve these unusual effects by means of tailored geo¬ 
metrical structures without relying on the intrinsic prop¬ 
erties of materials used in metamaterial fabrication. Con¬ 
ducting non-magnetic split-ring resonators (SRRs), and 
their derivatives, have been employed to demonstrate 
negative permeability by strongly coupling to the mag¬ 
netic field component of incident electromagnetic waves 
at microwave [5] or higher frequency bands PIT!. How¬ 
ever, when lower operational frequencies are needed, 
these resonator-based metamaterials require impracti- 
cally large dimensions to be able to compensate for 
ohmic losses [H]. One way to address this miniaturiza¬ 
tion problem is to use superconductors since they show 
very low surface resistance below their transition temper¬ 
ature mm- Moreover, tunability of the resonant fea¬ 
tures can be precisely achieved without lumped elements 
by means of temperature and magnetic field [TTHT5] . For 
instance, the plasma frequency of superconducting spiral 
meta-atoms have been tuned to the MHz range by means 
of temperature variation, and the plasmonic properties 
have been studied at essentially fixed frequency El- 

Of particular interest are radio frequency (rf) metama¬ 
terials where the free space wavelength Aq can be tens of 
meters. Despite the challenge of scaling the dimensions 
down, rf metamaterials offer a wide spectrum of applica¬ 
tions including high resolution magnetic resonance imag¬ 
ing [HHiT], magneto-inductive lenses [18] and wireless 
communication employing electrically small an¬ 


tennas. Antennas integrated with compact rf metamate¬ 
rials were speculated to have enhanced radiation power, 
high directivity and reduced return loss due to the unique 
electromagnetic properties of metamaterials [ 531121 ] ■ Zi- 
olkowski [^ [531 has proposed using a metamaterial with 
negative refractive index to boost the real radiated power 
of an electrically short dipole antenna (size Ao/1000) by 
almost two orders of magnitude compared to the dipole- 
in-free-space case. A large imaginary part (reactance) of 
radiation impedance Zrad=Rrad+iXrod of an electrically 
small antenna limits the radiation efficiency because free 
space has a purely real impedance |57|. When the an¬ 
tenna is placed in near-held proximity to a negative re¬ 
fractive index medium, the ^rad of the antenna can be 
cancelled out. Therefore the metamaterial acts as a nat¬ 
ural impedance matching network for the dipole, which 
maximizes the transfer of power from the antenna to free 
space in the far-held. The bandwidth and the operat¬ 
ing frequency of the antennas can also be reconhgured 
with tunable metamaterials, which is crucial for modern 
wireless communication systems relying on multi-band 
operation [53] . 

In this work, we demonstrate a tunable 3D rf meta¬ 
material composed of electrically small (^ Ao/2300) su¬ 
perconducting thin him resonators. The single resonant 
element is a planar spiral constructed by w= 3 /rm wide 
semi-circular strips with increasing radii. The center to 
center separation between two adjacent turns is s= 6 
/rm. The spiral is densely packed with a large number of 
turns, N= 200, where the inner and outer diameters are 
Di= 4.606 mm and Do= 7 mm respectively. The electro¬ 
magnetic response of the resonant spirals can be studied 
with a simple circuit model [28] ; i.e. an individual planar 
spiral can be treated as a magnetic resonator with an an¬ 
gular resonant frequency of a;o=l/vTC' where L and C 
are the total inductance and capacitance, respectively. 
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Possessing large geometrical inductance and capaci¬ 
tance (due to large number of turns) over a compact ge¬ 
ometry puts the individual spirals into a significantly sub¬ 
wavelength self-resonance regime. Therefore the meta¬ 
material formed by the 3D array of these spirals can be 
considered as a continuous effective medium character¬ 
ized with effective constitutive parameters; i.e. frequency 
dependent permeability /J,cfi(w) and permittivity eeff(a;). 
Previous theoretical [H] and experimental work 

on metamaterials with periodic circular conducting el¬ 
ements have shown that the fundamental resonant fre¬ 
quency of the metamaterial can be drastically lower than 
that of the single constitutive element mainly due to the 
mutual inductance (magnetic coupling) between the res¬ 
onant elements. However, electrical (capacitive) coupling 
can also play a role in the collective electromagnetic re¬ 
sponse of the metamaterial, thus interaction between the 
individual elements can be complex depending on struc¬ 
tural details of the metamaterial such as separation be¬ 
tween the resonant elements along the substrate, inter¬ 
layer spacing, array size and alignment/arrangement of 
the elements along the stacking direction |351 [33] . 

Here we study a 3D superconducting metamaterial 
with sub-wavelength spiral elements operating in the rf 
regime. We find that the mutual interactions between the 
individual spirals contributes to the effective response of 
the macroscopic 3D metamaterial. Moreover, we observe 
that both the center frequency of the resonant bandwidth 
in the rf transmission and the negative /Xeff(a;) region is 
highly temperature tunable. 


EXPERIMENT 

The metamaterial is patterned out of a 200 nm thick 
Nb thin film sputtered on to a 350 /im thick quartz sub¬ 
strate. We define circular spiral-shaped resonators [5] 
into the Nb film by using standard photolithography and 
reactive ion etching using a mixture of CF 4 /O 2 (10 % 
O 2 ). After fabrication, the quartz substrates are diced 
into 0.9 X 0.9 cm^ and 2.6 x 5.2 cm^ pieces to have in¬ 
dividual spirals and 2D arrays of 3 x 6 Nb spirals. To 
be able to extend the structure into three dimensions, 
four substrates containing 2D arrays are stacked along 
the c-axis of the substrate by using thin sheets of Roha- 
cell foam as spacers. The planes of the spirals are aligned 
in parallel and each substrate is separated from the next 
by 0.5 mm along the z-direction to produce strong cou¬ 
pling between the Nb spirals in the third dimension. The 
experiments have been conducted in an evacuated probe 
inside a cryogenic dewar and the temperature has been 
precisely adjusted by a temperature controller between 

4.2 K up to well above the superconducting transition 
temperature, ~ 9.2 K of Nb. 


RESULTS AND DISCUSSION 

Figurej^shows the rf-transmission |S'i 2 | through a sin¬ 
gle Nb spiral along with the 2D and 3D arrays of the 
coupled spirals. The measurements have been performed 
by means of two loop antennas sandwiching the samples 
(see the sketch in the inset of Fig. [^illustrating the 3D 
array of spirals between the antennas; the quartz sub¬ 
strates are not shown). The antennas are made of semi¬ 
rigid coax cables by shaping the inner conductors into 

1.2 cm-diameter loops and attaching them to the outer 
conductors. The other ends of the coaxial cables are con¬ 
nected to a 2-port vector network analyzer to measure 
the scattering matrix. The magnetic excitation of the 
sample is achieved via inductive coupling into the top 
(drive) loop and the transmitted signal is picked up by 
the bottom (pickup) loop. 

The loop antennas are 27.2 mm apart to minimize the 
direct magnetic coupling between them [Mj that can ob¬ 
scure the response of the 3D metamaterial. The single 
spiral demonstrates a sharp fundamental resonance peak, 
/o=a;o/27r at ^ 24.95 MHz and a dip at ~ 25.04 MHz at 

4.3 K with a loaded quality factor of 1273. The peak and 
the dip in the transmission spectrum correspond to con¬ 
structive and destructive interferences between the direct 
coupling of the loops and magnetic coupling of the loops 
to the metamaterial respectively |35|. We have observed 
higher order modes corresponding to shorter-wavelength 
rf-current standing waves in the spiral, which have been 
extensively studied elsewhere [HI [33H33| . Above the Tc of 
Nb, superconductivity is lost and the resonance is wiped 
out at 9.25 K. This demonstrates that such a compact 
device is not functional if made with thin films of nor¬ 
mal metals due to the large ohmic losses. An analytical 
model adapted to a similar experimental configuration 
showed that the transmission spectrum in the absence of 
the spiral resonator between the loops captures the ex¬ 
perimental transmission through the spiral taken above 
the Tc of the superconductor [35| . 

The low temperature resonance peak shifts to a 
significantly smaller frequency, 19.15 MHz in the rf- 
transmission spectrum of the 2D array of 3 x 6 spirals 
due to a significant change in effective parameters; espe¬ 
cially /ieff(a;) (33] . The small separation distance between 
the spirals ^ Ao/1875 strengthens the magnetic coupling 
between the individual spirals that leads to splitting of 
the resonant frequencies [30] . 

For 3D metamaterials with strong magnetic coupling 
in the stacking direction, the enhancement in capacitance 
due to interlayer interaction m further reduces the res¬ 
onant frequency at 4.3 K down to 18.83 MHz as seen 
in the transmission spectrum of the 3D metamaterial. 
The resonant features observed in the low temperature 
spectra of both 2D and 3D arrays of spirals disappear at 
temperatures above 9.25 K as expected. 
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FIG. 1: (Color online) Comparison of transmission |Si 2 | 
through a single spiral (green curve), 2D array of 3 x 6 spirals 
(red curve) and 3D array of the spirals formed by stacking 
four substrates including 2D arrays of spirals (blue curve). 
The lowest resonant frequency band decreases and develops 
fine splitting by adding more elements to the metamaterial. 
The maroon curve shows the data taken at the of Nb, 9.25 
K, for a single spiral, where the resonant peak disappears due 
to enhanced ohmic losses. 2D and 3D arrays also show the 
same feature when the sample is heated above 9.25 K. The 
inset shows a schematic of the 3D metamaterial sandwiched 
between two rf-loop antennas which are 27.2 mm apart from 
each other along the z-axis. 


By considering the effective medium approximation, 
the frequency dependent complex effective permeability 
of an array of non-magnetic conducting SRRs has been 
characterized with a Lorentzian function as [1] 


Meff(w) = Meff(‘^) + Veff(w) 


+ iToj 

where and are the real and imaginary parts 

of the effective permeability, F is the filling fraction, ojr= 
2TTfr is the angular center frequency of the resonance, and 
r is the damping term characterizing ohmic and radiation 
losses. The center frequency for the resonant bandwidth 
fr is determined by the resonant parameters of individual 
elements /o, as well as the mutual interactions between 
them [2^14^ . 

The wavelength of the electromagnetic excitation we 
have used is much larger than the dimensions (Dg, w, s) 
and separation of the spirals forming our metamaterial. 
Therefore, we can apply the effective medium approxima¬ 
tion and characterize the metamaterial with /ieff(a;) and 
eeff(u;). Since inductive coupling is quite prominent here, 
the former is the dominant effective parameter. Sharp 
dips of the transmission spectra just above the resonant 
frequencies in Fig. (below the background signal set 


by the normal state data taken at 9.25 K) imply the 
presence of negative permeability. To be able to char¬ 
acterize the negative magnetic response of our metama¬ 
terials, we have retrieved from the experimental 

scattering matrix m collected by two loops connected 
to a network analyzer (see the inset of Fig. 1). The in¬ 
duced electromotive force in the pickup loop as a func¬ 
tion of frequency due to the magnetic coupling of top 
loop can be written as emf(a;)ref= -iujA^qHz when there 
is no metamaterial between two loops and emf(w)nieta= 
when there is metamaterial sandwiched 
between them. Here, A is the area of the pickup coil, fj,o 
is the free space permeability and Hj, is perpendicular 
component of the magnetic field. By combining these 
two equations one can approximate the complex effective 
permeability as a ratio, /ieff(a;)=emf(w)meta/emf(a;)ref 
<S' 2 i(w)meta/<S' 2 i(w)ref- Here it is assumed that 77^, is the 
same in both measurements [41]. 

Figure]^ shows the real (red curves) and imag¬ 

inary /i"g(a;)(blue curves) parts of the retrieved perme¬ 
ability for (a) the single spiral and (b) the 3D metamate¬ 
rial composed of similar spirals. The overlapping dashed 
curves are the Lorentzian fits according to Eq. (1) which 
reasonably capture the retrieved data at 4.3 K. The fit 
parameters are uJoj2n= 24.955 MHz, T/2'k= 0.03 MHz 
and F= 0.0068 for the single spiral and a;j./27r= 18.828 
MHz, F/27r= 0.064 MHz and F= 0.012 for the 3D array. 
The enhanced loss in the 3D array compared to an in¬ 
dividual spiral can be explained in part by the stronger 
inductive coupling to the excitation loops |43| and the 
excitation of magneto-inductive waves in the metamate¬ 
rial |31|. We observe that the smooth well-defined reso¬ 
nance of the single spiral splits into multiple resonances 
with magnetic coupling into other spirals in the 3D array 
(see the insets of the plots showing the resonant features 
in detail). Just above the resonant frequency, the 
becomes negative in a narrow frequency band. These for¬ 
bidden bands of frequencies range from 24.958 to 25.038 
MHz for the single spiral and from 18.838 to 18.924 MHz 
for the 3D metamaterial at 4.3 K. 

Having discussed the low temperature magnetic re¬ 
sponse of the metamaterial, we now turn to the tunability 
of fj,eg{uj) as the metamaterial is heated up to a temper¬ 
ature above the Tc of Nb. Figure |^a) and (b) show 
the temperature dependence of the real and imaginary 
parts of the retrieved /ieff(w) for the 3D metamaterial. 
The narrow frequency band where the is negative 

shifts to lower frequency as the temperature increases as 
seen in Fig. da). The same effect is also seen in the 
peaks of the that are often associated with the 

resonant frequencies [Til liS] (see Fig. ib)). The fre¬ 
quency shift in the spectra is a result of a decrease in the 
density of superconducting electrons (superfluid density) 
with temperature UsiT) which is inversely proportional 
to the square of magnetic penetration depth Al(T). As 
the fields of the electromagnetic waves penetrate into the 
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FIG. 2: (Color online) Real and imaginary 

parts of the effective complex permeability retrieved from the 
experimental scattering matrix at 4.3 K for (a) a single spiral, 
(b) 3D array of such spirals. The dashed maroon and cyan 
curves are Lorentzian fits to the retrieved data according to 
Eq. (1) demonstrating a significant agreement. The insets are 
shorter frequency range plots of the data showing the details 
in the vicinity of the resonant frequency. A single smooth 
resonance peak of the individual spiral splits into multiple 
resonances for the 3D metamaterial. 


metametarial deeper with increasing temperature, the to¬ 
tal inductance of the individual spirals increases due to 
the enhanced geometrical and kinetic inductance, Lk oc 
Xl{T) coth[t/AL(T)] where t is the thickness of the Nb 
film used to make the metamaterial [46]. The enhanced 
total inductance leads to a reduction in the resonant fre¬ 
quency m- Moreover, transmission decreases in mag¬ 
nitude and is smeared as is approached due to the 
enhanced ohmic losses arising from normal carriers, char¬ 
acterized by r in Eq. (1). 

The narrow frequency band in the real part of /ieff(w) 
possesses negative values up to ~ 8.3 K beyond which it 
starts to show near zero and then positive values. This 
crossover temperature corresponds to a similar value at 
which we start to see a discernible downward trend in 
the temperature evolution of fr{T) and an upturn in 
the damping term r(T). To study this behavior, we ex¬ 
tracted ojr (thus fr) and T as a function of temperature 
through fitting the retrieved /iefi(a;) of the 3D metama¬ 
terial according to Eq. (1). As seen in the main panel of 
Fig. El for a long range of temperatures below T^, both 



Frequency (MHz) 

FIG. 3: (Golor online) Temperature evolution of (a) real 
/igff(tj) and (b) imaginary /i"fj(a;) parts of the retrieved per¬ 
meability for the 3D metamaterial. Both the frequency band 
where the fJ^giuS) is negative and fr shift to smaller values of 
frequency with increasing temperature due to the enhanced 
inductance and losses. 

fr and r are quite constant since there is no significant 
change in ns{T) E5]. As temperature increases, ns(T) 
starts to decrease leading to a weaker response in oppos¬ 
ing the external magnetic held. Beyond 8.3 K, we see a 
slight decrease (increase) in fr (T) until the temperature 
is very near Tc- In the vicinity of Tc, we observe a steep 
drop (increase) in fr (F) mimicking the behavior of rigifT) 
{nN(T), normal fluid density) in this temperature regime. 
Here, Ai,(T) and thus the total inductance dramatically 
increases (i.e. Lk diverges) and accordingly ns{T) ap¬ 
proaches zero. At Tc, ohmic losses are enhanced due to 
the dominance of n]\[{T) over ns(T) and resonant charac¬ 
teristics are completely wiped out. To demonstrate this 
more clearly and to test the validity of our extraction, we 
show l-[friT)/fr{Tmin)]'^ VS. r(r)/r(r„^i„) in the inset 
of Fig. El Both quantities are roughly proportional to 
the normal huid density n^(T) so the plot should be ap¬ 
proximately a straight line. Though we see a nonlinearity 
in the vicinity of the base temperature (due to the near 
absence of a frequency shift), the trend is quite linear at 
the higher temperatures, as expected. 

SUMMARY 

We have developed a low-loss and ultra sub-wavelength 
3D superconducting metamaterial. The engineered com- 
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FIG. 4: (Color online) Temperature evolution of the extracted 
resonance frequency fr (blue squares) and damping factor F 
(red diamonds) by fitting the retrieved permeability of 3D 
metamaterial according to Eq. (1). In the vicinity of Tc, fr 
shows a steep drop as F dramatically increases. The inset is 
a plot showing the quantities expected to be proportional to 
the normal fluid density UNiT). 

pact metamaterial shows a strong magnetic response to 
electromagnetic waves. The resonant frequency and fre¬ 
quency band with negative /ieff(a;) are precisely tuned by 
means of temperature without lumped elements. We ob¬ 
serve a crossover from negative to positive at a 

temperature near T), at which a dramatic increase occurs 
in losses (F) signaling a significant decrease in superfluid 
density ns(T). The work contributes to ongoing research 
on tunable negative permeability and near zero perme¬ 
ability metamaterials which can play an important role 
to improve the performance and directivity of ultra-small 
rf antennas. 
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